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ABSTRACT 

Based on the largest homogeneously reduced set of Type la supernova luminosity data 
currently available - the Union2 sample - we reconstruct the expansion history of the 
Universe in a model-independent approach. Our method tests the geometry of the 
Universe directly without reverting to any assumptions made on its energy content 
and thus allows us to constrain Dark Energy models in a straightforward way. This is 
demonstrated by confronting the expansion history reconstructed from the supernova 
data to predictions of several Dark Energy models in the framework of the tuCDM 
paradigm. In addition, we test various non-standard cosmologies such as brancworlds, 
f(R) and kinematical models. This is mainly intended to demonstrate the power of the 
method. Although statistical rigor is not the aim of our current study, some extreme 
cosmologies clearly disagree with the reconstructed expansion history. We note that 
the applicability of the presented method is not restricted to testing cosmological 
models. It can be a valuable tool for pointing out systematic errors hidden in the 
supernova data and planning future Type la supernova cosmology campaigns. 

Key words: supernovae: general - cosmology: cosmological parameters - cosmology: 
observations - cosmology: theory 
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Type la supernovae (SNe la) are currently t he best (relative 
distance indicators out to redshifts of z ~ 1 (Tammannll97T 
Colgate!! 19791: iRiess. Press fc Kirshnerlll99a ISchmidt et al .1 



1998; Pcrlmuttcr ct al. 



1999! ). 



As such, they have been 
instru mental in shapin g our current picture of the Uni- 
verse |LeibundgutJl200ll . I 20081 '). In particular, the notion of 
an accelerated expansion rate of the Universe was estab- 
lished a decade ago based on SN la dista nce measurements 
(|Riess et al.ll 19981 ; IPerlmutter et alj|l999h . 

For lack of a deeper understanding, the cause of this 
acceleration is commonly parametrized in standard ACDM 
cosmology as a Dark Energy component that currently 
dominates the energy contents of the Universe (see, e.g., 
I Turner fc Hutererll2007l . for a recent review) . The most plau- 
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sible explanation for it may be the vacuum energy (or 
cosmological co nstant) with a constan t equation of state 
w = p/p = -1 (jTurner fc White! Il997h . In this case, Dark 
Energy would be an elastic and smooth fluid exerting a 
repulsive g ravity that produces the o b served accelerated 
expansion (lOstriker fc Steinhardtlll995l ; iLiddle et al]|l996l ; 
iTurner fc Whiteill997h ~ 

So far, however, all attempts to compute the vac- 
uum energy have led to values that are abo ut 55-120 
order s of magnitu de off the observed value (|Weinberd 
1 19891 ; ISahnil 120021 '). Other possibilities allow for a Dark 
Energy equation of state varying with time and avoid 
the cosmic coincidenc e and the previously mentioned 
fine-tuning problems (IZlatev. Wang fc Steinhardtl 1 19991 : 
iFrieman. Turner fc Hutererll200c€ 



These concepts describe the vacuum energy to be a dy- 
namical, evolving scalar-field slowly r olling toward its low- 
est energy state (|Freese et alj 1 19871 ; IFrieman et al.l 1 19951 : 
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Wetterichlll995l;lTurner fc Whitdfl997l ; ICaldwell et ai1ll99Sl ; 



Steinhardt 19991 ). In this context, the simplest parametriza- 



tion of the Dark Energy equatio n of state w(a) depending on 
the cosmic s cale factor a reads (|Chevallier fc PolarskiluOOll ; 
lLinder!l2003l ): 



w(a) — wo + ra a (l — a), 



(1) 



where wo is the present-day value of the equation of state 
and w a accounts for its time-dependence. 

Alternative approaches describe cosmic accelera- 
tion as a manifestation of new gravit ational physics 



rather than Dark Energy (see, e.g ., Deffavet) 200 ll; 



Deffavet et al"1l2002l; ICarroll et aLlfeoO^ ; iNoiiri fc Odintsovl 
20061 : lAmendola et al.l 1200 7l ). Instead of adding an extra 
term to the energy-momentum tensor on the right side of 
Einstein's equations, these modify the geometry terms on 
their left side in order to reproduce the observations. 

Hence, we face the situation that apart from ACDM 
a large variety of cosmological models has been proposed 
to account for cosmic acceleration. Testing their validity is 
an important but challenging task. The present-day nearby 
Universe is fairly well-known and therefore all models under 
consideration reproduce its characteristics, or contain free 
parameters that can be tuned to do so. This degeneracy is 
difficult to break with currently available data. 

Usually, Dark Energy models are constrained by start- 
ing out from a Friedman cosmology in which the expansion 
function H is parametrized in terms of the contributions of 
radiation (r), matter (m), curvature (fc), and Dark Energy 
(de) to the energy density as 



iDalv fc Diorgovskil (|2003l . 120041 ) combining SNe la luminos- 
ity dista nces with angular - diameter di stances from radio 
galaxies. ISeikel fc Schwarj i|200Sl . |2009| ) tested the signifi- 
cance of cosmic expansion directly from SN la data in a 
model-independent way. 

The goal of this work is to apply the method of 
iMignone fc Bartelmannl (|2008l ) to the most complete SN la 
data set currently available and to constrain some flavors of 
Dark Energy models. This is intended as a demonstration of 
the general power of SN la distance measurements for testing 
Dark Energy models when analyzed in a model-independent 
way. A rigorous statistical treatment that would let us rule 
out specific models is beyond the scope of this paper. 

The paper is organized as follows. In Section [2] the es- 
sential aspects of the model-independent methodology are 
reviewed. The application of the method to luminosity- 
distance measurements is discussed in Section [3] A com- 
parison between SN la data and several Dark Energy mod- 
els is presented in Section [4] In Section [5] the predictions 
for the expansion history of the Universe of several non- 
standard cosmologies are confronted with our reconstruc- 
tion from SN la data. Ways of improving the reconstruction 
with new data that may become available through current 
surveys are pointed out in Section This shows that our 
model-independent analysis can be instrumental in planning 
observational campaigns. Moreover, it can also point to po- 
tential systematics in the data that do not generally arise 
from traditional ways of analysis. This aspect is highlighted 
in Section [7] Finally, conclusions are drawn and future per- 
spectives are discussed. 



H 2 {a) 



= H 2 E 2 (a) 



= m 



—r H 5 9~ + "deO-r (a) 

a 4 a s a 2 



(2) 



where Ho and f2...o denote the Hubble constant and the 
present-day density parameters corresponding to the differ- 
ent components, respectively. A possible time-dependence of 
the Dark Energy equation of state is captured by the func- 
tion F(a) 



(3) 



F(a) = exp{-3^ Q [l + ™(a')]^}. 

The cosmological parameters occurring in Eq. ([2]) are de- 
termined from fits to observations and Dark Energy models 
are usually assessed by confronting their predictions to these 
parameters. The significance of this approach is, however, 
limited since it automatically assumes a Friedman model. 
When working with SNe la this constraint is unnecessary as 
they probe the geometry of the Universe directly and no as- 
sumptions on the form of the energy-momentum tensor are 
required to derive the expansion history of the Universe. 

Con sequently, our analysis follows a recently developed 
method (|Mignone fc Bartelmannl |2008l ) to reconstruct the 
expansion history of the Universe in a model-independent 
fashion from luminosity distance data. The idea of a model- 
independent reconstructio n extracted st r aight f rom the data 
was already proposed in IStarobinskvl (Il998|) and recon 
structi ons of this kind were carried out by 
IShafielool (|2007l ) using SN la data, by 



S hafieloo et al 



Fay fc Tav akol 



200' 

20061 ) adding constraints from measurements of baryon 
acoustic oscillations (BAO) to the SN la data, and by 



2 MODEL-INDEPENDENT METHOD 

The minimal a ssump tions adopted by 

IMignone fc Bartelmannl (|2008l ) are that the expansion 
rate is a reasonably smooth function and that the Uni- 
verse is topologically simply connected, homogeneous and 
isotropic, i.e. it is characterized by a Robertson- Walker 
metric: 



ds 2 = cdt - a 2 (t) [d X 2 + fl{x){dQ 2 + sin 2 8d<j> 2 )] 



(4) 



The Robertson- Walker metric allows us to define an 
angular-diameter distance by 

Da (a) = af k \x(a)l (5) 
with the comoving angular-diameter distance 

(sinx (k = 1, spherical) 
X (k = 0, flat) (6) 
sinhx (fc = — 1, hyperbolic) 

and the comoving distance 

nl dx 



(7) 



H J a x*E(x) 

Through Etherington's relation dEtheringtonl 1 19331 ') . 
which holds for any space-time, we can relate the luminosity 
distance to the angular- diameter distance: 



D L (a) = -^D A {a) 



(8) 



This allows us to write the former as an integral of the in- 
verse of the expansion rate 
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D L (a) = -^i 
Ho a 



d.r 



c 1 
H a 



1 dx 



(9) 



x 2 E{x) 

with e(a) = E~ (a) and fc = set in Eq. ((6]) for simplicity 
of notation. It is worth noticing that the choice of fc = 
does not affect the fundamental method and can be dropped 
without change of principle if needed. The Hubble length 
c/Ho, by which the luminosity distance is scaled, shall be 
dropped in the following for the sake of brevity. 

From Eq. (|9]) the derivative with respect to a is taken 

l %e{x)- e -^. (10) 



DUa) 



This expression can be transformed to a Volterra integral 
equation of the second kind for the unknown function e(a), 



e(a) = —a 3 D' L (a) + A / ^e(x 



(11) 



with the inhomogeneity f(a) = — a 3 D' L (a) and the simple 
kernel K{a,x) = x~ 2 . The general parameter A will later 
be fixed to A = a. For now it is introduced to make the 
connection to a class of equations for which solutions are 
known to exist and to be uniquely describ ed in terms of a 
Neumann series (see lArfken fe Weberill995h : 



e(et) = y~] X l ei(a). 

i=0 

A possible choice for the function e, would be 

eo(a) = /(a), 

e n (a) = / K(a,t)e n -i(t)dt, 



(12) 

(13) 
(14) 



where, for the guess of eo, the approximation of either the 
integral or A to be small has been made in Eq.{TT]). This ap- 
proximation is valid in all important cosmological cases and 
it is subsequently improved until convergence is reached. 
It essentially means that starting, for instance, from the 
ACDM cosmology observations say that deviations must be 
small, if they exist at all. 

Equation ()lip involves the derivative of the luminosity 
distance with respect to the scale factor a. Observations of 
SNe la provide measurements of the distance modulus, fn, 
and redshifts, 2, (or scale factors at = (1 + Zi)^ 1 ). Thus 
one could think of taking the derivative of the luminos- 
ity distance directly from the data. However, this is an in- 
convenient procedure since the result would be extremely 
noisy and the determination of D[ j (a) would be unreliable. 
Therefore, it is necessary to suitably smooth the data in 
the first place by fitting an adequate function Dl(£s) to 
the measurements D\ J {ai). The derivative in Eq. (|11[) can 
then be approximated by the derivative of D^(a). Thus, the 
derivative of the fitted data is taken as a representation of 
the derivative of the real dat a. The method proposed by 
iMignone fc Bartelmannl (|2008T ) achieves this goal via the ex- 
pansion of the luminosity distance Dl {a) into series of (in 
principle) arbitrarily chosen orthonormal functions Pj(a): 

j 

£>L(a)=X)ctfy(a). (15) 
j=o 

Suitable orthonormal function sets can be constructed by 
Gram-Schmidt orthonormalisation from any linearly inde- 
pendent function set. The J coefficients Cj are those which 



minimize the \ 2 statistic function when fitting to the data. 
Therefore, with this representation of the data, the deriva- 
tive of the luminosity distance function is simply given by 

j 

3=0 

Due to the linearity of Eq. (|lll) . it is possible to solve it 
for each mode j of the orthonormal function set separately. 
Thus, the final solution in terms of a Neumann series is 



e(o) 



(j) 



(a), 



(17) 



that is, the measured coefficients of the series expansion give 
the solution for the expansion function. 

The number of terms to be included in order to pro- 
duce an acceptable fit to the data depends on the choice of 
the orthonormal basis. Due to the limited quality of current 
observational data, however, it is expected that the coef- 
ficients of this expansion can only be determined up to a 
certain order. Therefore, although the basis of the expan- 
sion is arbitrary with ideal data, it is not in practice. A 
preferred choice would minim i ze the number of required co- 
efficients. iMaturi fe Mignond (|2009r ) suggested an optimal 
basis system derived from principal component decompo- 
sition of cosmological observables. The use of this optimal 
basis would certainly reduce the coefficients needed to ob- 
tain an accurate reconstruction and we shall do so in future 
work. In the current analysis, howe ver, we settle for the ar- 
bitrar ily chosen basis proposed by IMignone fc Bartelmannl 
which uses the linearly independent set 



Uj(x) 



(18) 



orthonormalized via the Gram-Schmidt method. The or- 
thonormalization interval is chosen to be [a m i n , 1], where 
a m in = (1 + Zmai) -1 is the scale factor corresponding to 
the maximum redshift z ma x in the supernova sample. This 
way, an arbitrary set of orthonormal functions Pj (a) is con- 
structed. When projecting the luminosity distance Dl(o) 
onto these basis functions, we can solve for the expansion 
coefficients. In this case, at least the first five coefficient s 
are different from zero (see IMignone fc Bartelmannl 120081 '). 
However, as shown in Table [1] only the first three can be 
determined with current data; the fourth and fifth coeffi- 
cients lose significance. With future space-based telescopes, 
new and better-quality data will be become available for this 
kind of analysis. This should procure an improvement on the 
accuracy of the reconstruction by allowing for determining 
more coefficients in the expansion. 



3 APPLICATION TO SN IA DATA: THE 
UNION2 SAMPLE 

We apply the model-independent method to the largest 
homogeneously reduced SN la sample currently avail- 
able. The Union2 sample (|Amanullah et ail l2010f ) con- 
sists of 557 SNe la. It includes the recently extended 
dataset of distant supernovae observed with the HST 
|Riess et all 120071; [A manulla h et all l20ld '). the data from 



the SNLS (lAstier et all 



l2007i ; IWood-Vasev et al 



2006), ESSENCE (iMiknaitis et al 



2003) and SDSS (|Holtzman et al 
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20081) surveys, sev eral compilations from literature (e.g. 



Hamuv et all 1 19961), and the new data from nearby SNe la 
of lHicken et all (|2009l ). 

In Fig. [T] the Union2 sample is shown, together with 
the best fit of the luminosity distances to Eq. (|15[) when 
using the first three terms in the expansion (dotted line, see 
Table[T]for values). There is considerable scatter around the 
fit, mainly introduced by the distant SNe which are more 
difficult to calibrate and are usually afflicted by host-galaxy 
extinction and survey selection effects. 

In the following two sections, we compare the expansion 
history derived from the Union2 sample with predictions of 
various cosmological models. This comparison, however, is 
not fully consistent. SNe la are not standard candles but 
have to be calibrated as distance indicators applying an em- 
pirical calibration. This calibration procedure introduces a 
dependence between the calibrated measurements that gives 
rise to finite non-diagonal entries in the covariance matrix 
for the distance modulus. 

Moreover, the parameters of the calibration are deter- 
mined simultaneously with the cosmological par ameters in 
the SALT2 light curver fitter (|Guv et all l2007t ) assuming 
a ACDM model. As the Union2 sample was calibrated with 
SALT2, this applies to the data used in the following and im- 
plies that the parameters are model-dependent. This clearly 
contradicts the intended model-independence of our method 
and would prevent precise fits to non-ACDM cosmologies. 
Such fits, however, are not our goal here. Apart from that, 
ACDM happens to be an extraordinary good fit to the data 
and models with more degrees of freedom do not improve 
the quality of the fit more than statistically expected. That 
means that the correction parameters will not change signif- 
icantly if one assumes a different model, as long as it roughly 
matches ACDM. Models that show a H(z) drastically dif- 
ferent from ACDM could also have different parameters in 
the calibration, but are not interesting as they would not 
agree with the data. Additionally, the calibration parame- 
ters should be redshift independent and therefore, essentially 
decoupled from cosmology. Consequently, we argue that the 
induced covariance is small and the variation of the param- 
eters with cosmology is negligible for our purposes. 



4 COMPARISON TO ACDM AND DARK 
ENERGY MODELS 

4.1 ACDM cosmology 

The result of our reconstruction of the expansion rate H (a) 
from the Union2 sample is compared with a ACDM model 
in Fig. [2] Here, we adopt the values for the today's density 
parameters from t he bes t fits to the Union2 sample given by 
lAmanullah et al] (|2010h - Q m = 0.27A±^° t°. The Hubble 
constant is assumed to be h = 0.7 (Ho = lOO/i km s" 1 
Mpc -1 ) and the equation of state to be constant (w — —1). 

Thanks to the improved quality and size of the Union2 
data se t, the reconstruction ha s smal ler errors bars than 
that of iMignone fc Bartelmannl (|2008t ). who made use of 
the SNLS data only. The ACDM model is in good agree- 
ment with the supernova data, although its slope is slightly 
different. This leads to a deviation at intermediate values 
of a, but, within the error bars, the ACDM model is still 
consistent with our reconstruction. 



Union2 Sample 



SNLS 
ESSENCE 

Samples from literature 
New Nearby - Hicken 2009 
HST - Riess 2007 
New HST - Amanullah 2010 
SDSS - Holtzman 2008 
Fit for D L 




Figure 1. Luminosity distance plotted against scale factor a for 
the the Union2 sample. The dotted line is the best fit for the 
luminosity distance obtained by fitting the data to Eq. 11151 and 
using three expansion coefficients. 



Union2 Sample 



3a reconstruction 

ACDM model: 

fl m = 0.274 + ° D 3 4 ° (best tit from Amanullah et al. 201 0) 




Figure 2. The reconstructed expansion rate, with 3<r errors, ex- 
tracted from the Union2 sample. The red lines represent a ACDM 
model assuming the bes t-fit cosmological parameters reported by 
lAmanullah et al] l |2O10h . It is shown here only for the sake of 
comparison, this model was not used in our calculations. The 
bottom panel shows the residuals between the reconstruction and 
the model. The dotted lines are the limits of the model varying 
Q m with H fixed to 70 km s -1 Mpc -1 . 



4.2 Dark Energy models 

Relaxing the assumption of constant equation of 
stat e, we analyze a Dar k Ene rgy model proposed 
by iRapetti. Allen fc Welieil (|2005h . It is an exten- 
sion of the parame t rizatio n in Eq. CO prop osed by 
IChevallier fc Polarskil (|200lj ) and iLindeiT (|2005f l. which 
assumes a fixed transition redshift (z t = 1) between the 



current value of the equation of state and the value at early 
tim es, w et = wp + w\. In contrast, the model discussed 
by IRapetti et al.1 (l200fih introduces z t as an extra free 
parameter so that the equation of state w can be written as 
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Table 1. The first five coefficients of the reconstructed expansion. 
Coefficients Fitted Values 



Union2 Sample Vs. DE Models 



CO 
CI 
C2 
C3 
C4 



2772 ± 14 
-1389 ± 15 

73 ± 11 
12± 5 

-3 ± 1 



w(a) 



Wctz + w zt w ct (l - a)a t + w (l - a t )a 



Z + z t 



a(l - 2a t ) + a t 



(19) 



where at is the transition scale factor. From the Friedmann 
equation, the expansion function in terms of redshift is then 
given by 



H{z) = H v/fi m (l + 2)3 + fW(z) + n k (i + zY, (20) 
with 

f( z \ — Q _|_ z ^3(l + uj t) e -3(iu t-u)o) 9 (z;zt). 

Here, the function g(z;z t ) is defined as 



g{z;z t ) 



Zt 



hi 



Zt 



(22) 



Zt - l \z t + z i 

iRapetti et all (|2005l ) constrained the best-fit cosmolog- 
ical parameters for different possibilities within this model, 
varying the number of free parameters. The y made use of 
SN la (the Gold sample bv lRiess et al.ll2004h . X-ray galaxy 
clusters and CMB data for this analysis. In Fig. [3] we 
compare those possibilities to the model-independent re- 
construction of the expansion function extracted from the 
Union2 sample. This is meant to illustrate that standard cos- 
mologies within the toCDM paradigm with different choices 
of parameters do not fit the data in all cases. This is not easy 
to see in other analyses which are model-dependent. Our re- 
construction of H(z) can be a useful tool since it allows us 
to rule out cosmological models based entirely on the data. 
For example, the two most extreme models shown in Fig. [3]- 
one with a fixed z t = 0.35 (green dashed-dotted line) which 
splits the SN and the cluster data sets into similarly low 
and high redshift subsamples, and the one with arbitrary 
Zt (blue dashed line) - are clearly inconsistent with our re- 
construction. The other models (marginally) agree with the 
SN measurements within the error bars. However, they have 
slightly different slopes and - as for the ACDM model - some 
tension exists at intermediate and small values of a (inter- 
mediate and high redshifts). 



5 GOING BEYOND ACDM 

Despite the good agreement with cosmological observations, 
the ACDM model is conceptually problematic, as it leaves 
the fine tuning and coincidence problems unresolved. A 
number of alternatives to the ACDM parameterization have 
been suggested in order to address those difficulties. We 
study some of these alternatives here and compare them 
to our reconstruction. 



3c reconstruction 
DEI(Constanttlat): w n =-1.05 
DE2(flat, z,=1 ): w =-1 .1 0, w 1= 0.23 
DE3(split z,=0.35): w =-1 .30, w,=0.6 
DE4(arbitrary z,): w =-1 .27^=0.61 




Figure 3. The reconstructed expansion history, with 3a errors, 
extracted from Union2 sample and confronted with several Dark 
Energy mo dels with the best fi t cosmological parameters con- 
strained bv lRapetti et al" 1 1 |2005|) . 



5.1 Modified Gravity 

The idea of modified gravity differs from dark energy mod- 
els. It does not describe the late-time acceleration as being 
caused by some unknown energy component but suggests 
that General Relativity may be inaccurate at large scales. 
Such scenarios are likely to diverge in their predictions for 
H(z) from ACDM and therefore they can be constrained 
from the reconstructed expansion history of the Universe. 



5.1.1 The f(R) models 

One way of modifying gravity theory is by adding terms 
to the Ricci scalar (R) in the Einste in-Hilbert Lagrangian 
(|Starobinskvll 19801 ; ICarroll et al.|[200i) . The Palatini formal- 
ism provides second order differential field equations that 
can account for the pr esent cosmic acc eleration without any 
need of Dark Energy (|Fav et alj|2007l ). Several parameteri- 
sations of f(R) are possible. For the sake of simp licity, we 
follow the parametrization of lCarvalho et all (|200a) here (for 
detai l s on other para meterisations see, e.g., Hu fc Sawickil 
120071 ; iFav et all 120071 '). The general functional f(R) is as- 
sumed to have the form 

£_ 

R n ' 

where R is the Ricci scalar; n, j3, and Q m , are the parameters 
of the model. ACDM is recovered for n = The expansion 
rate can be written as 



f(R) = R 



(23) 



H 2 {z) = Hi 



m m0 (i + z) 3 + f/Hi 



6/'C 2 



with 



9 /" Hin m0 (i + z) 



(24) 



(25) 



2 /' Rf" - f 

where the notation /' = df/dR, f" = d 2 f/dR 2 has been 
adopted. By evaluating these equations at z = it is possible 
to obtain a relationship between n, /3, and Sl m . Thus, if we 
specify the values of two of them the other is immediately 
fixed. In our analysis we treat n and Q m as free parameters. 
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Union2 Sample Vs. f(R) Models 



Friedman equation is a modification of the Generai Rela- 
tivistic case and reads 











I 


3a reconstruction 






n=0.027, p=4.63, il m =0.28 (Amarzguioui et al. 2006): — 




— 


n=-0.09, p=3.6, fl m =0.27 (Fay et al. 2007): 






n=0.03, 0=4.7, a m =0.26 (Carvalho et al. 2008) — 




1 





Figure 4. The reconstructed expansion history, with 3cr errors, 
extracted from Union2 sample compared with several f(R) mod- 
els. The values for the models are taken from previously reported 
fits to BAO and CMB data in combination with independent 
H(z) (black dashed-dottcd line) and luminosity distance mea- 
surements (red solid — SNLS sample - and blue dotted — Gold 
sample - lines). 



We compare the model-independent expansion history 
derived from the Union2 sample with different f(R) mod- 
els in Fig. [4] The adopted values for the model parame- 
ters are taken from pr eviously reported fits in the literature. 
ICarvalho et al.l f|200Sh obtained constraints on the (n, fi m ) 
parameter space from the combination of BAO and CMB 
data with independent determinations of the Hubble param- 
eter at different redshifts. This model (black dashed-dotted 
line) has a substantially different slope and diverges from 
our reconstruction, 

In similar studies, other authors have combined SNe 
data with BAO and CMB. We consider here the number s 
obtained bv lFav et all |2007l ) and lAmarzguioui et all (|2006l ) 
using the SNLS and Gold samples, respectively. We find 
these models (red solid and blue dotted lines) closer to our 
reconstruction although they present deviations at high red- 
shift (for the Gold sample) and at intermediate redshifts (for 
the SNLS sample). 



5. 1.2 Brane worlds 

Another possibility to modify gravity is to consider higher 
dimensional models. Braneworld cosmology describes the 
observable Universe as a three-dimensional brane embedded 
in higher dimensional spacetime (or "bulk" ) . All matter par- 
ticles and fields are trapped on the brane (i.e. conservation 
of mass-energy holds firm) while gravity is free to propagate 
in the brane as well as in the bulk. 

There are di fferent scenario s within brane-world cos- 
mology (see, e.g.. lMaartensll2003 : lAlam fc Sahnill200l ). We 
consider here, f or example, the Dvali-Gabadadze- Porrati 
(DGP) model (|Dvali. Gabadadze fc Porratil |2000| ). This 
model allows the extra dimension to be large, and its gener- 
alisation to a Friedmann- Ro bertson- Walker brane produces 
a self-accelerating solution (|Deffavet|[200ll ). The resulting 



„a H 8irGp 
ti = — - — , 



where the cross-over scale r c is defined as 
1 



(26) 



(27) 



# (l-«m)' 

The additional term H/r c in Eq. (|26[) behaves like a Dark 
Energy component with an effective equation of state that 
evolves from w — —1/2 for z > 1 to u = - lin the distant 
future. 

For a toCDM model (assuming that Eq. |l} holds) it is 
possible to obtain an expression f or w(a) which mimics the 
evolu tion of the DGP model (see iTang. Weller fc Zablockil 
l2006h : 



w(a) 
to 



[(r CJ f7 )- 1 + 27 ? ]7 ? ' 



(28) 



va 



3 (o=l) = -: 



+ 1/2 (r c Ho) 3 . At present-day, this tends 



(29) 



■a. 



This modification allows for a description of the expan- 
sion history as well as of the growth of large-scale structure. 

One can even go beyond DGP as an isolated theory 
and consider a phenomenological model which is motivated 
by the concept of an extra dimensio n with infinite extent . 
This is the so called mDGP mo del (|Dvali fc Turner! 120031 ; 
iThomas. Abdalla fc Wellerl 120091 ') . It interpolates between 
ACDM and the DGP model and allows for the presence of an 
extra dimension through corrections to the Friedman equa- 
tion by introducing a parameter a. The modified Friedman 
equation reads 



H 1 



H° 



8nGp 



(30) 



with the cross-over scale r c defined as 

rc = (l-ftn)^H _1 . (31) 

In Fig. [5] the comparison of our model-independent 
reconstruction with some DGP and mDGP models is pre- 
sented. The values of the model parameters are taken from 
previously reported fits to actual cosmological data. 

We evaluate two pure DGP mo dels (a = 1) ado pting 
the best-fit values fo r f2 m reported in lGuo et al.l (|2006l ) and 
iLiang fc Zhul (|201ll ) (green dash-dotted and black dashed 
lines, respectively). In the first case, the constraints were 
obtained from the Gold and SNLS samples in combination 
with BAO. In the second, the best fit values are found by 
combining cosmology-independent Gamma Ray Burst and 
SNe la data, with BAO, CMB and H(z) measurements. The 
latter constitutes so far the strongest constraint obtained for 
the DGP model. 

W e also consider the results obtained bv lThomas et all 
|2009h adding weak-lensing data to BAO and SNe (red dot- 
ted line). In that particular study, the authors find an upper 
limit for the a parameter (a < 0.58 at 68 per cent confi- 
dence level), but are not able to give constraints on the n m 
parameter. Therefore, and to better understand the effect of 
changing f2 m for a given a, we make here use of the best-fit 
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Figure 5. The reconstructed expansion history, with 3<r er- 
rors, extracted from Union2 confronted with several DGP/mDGP 
models. The parameters of the model are taken from previous re- 
ported fits to several cosmological data. 



Figure 6. The reconstructed expansion history, with 3a errors, 
extracted from the Union2 set compared with several kinematic 
model s with best-fit parameters obtained by Guimarac s et all 

teoogh . 



values reported in Liang fc _Zhul l|201ll ) for the DGP model 
and in lAmanullah et all |2010l) for the ACDM model. 

For the sake of comparison, a ACDM model correspond- 
ing to a = (blue solid line) and a model with negative a 
(grey long-dashed line) are also included in Fig. [5] 

As it has been found in prior studies, (see, for exam ple, 
iMaartens fc Maierottol l200fj ; iFairbairn fe Goobarl l2006h . a 
pure DGP cosmology is disfavored by the SNe data. Mod- 
els with negative values of a also disagree with our recon- 
struction. However, it is still early to break the degeneracy 
betweeten ACDM and mDGP models with < a < 0.50. 



jerk parametrization, j(z) — jo- The last model, M4, derives 
from an expansion of the luminosity distance with free go 
and jo parameters resulting in 



Dl(z) 



c 

1 

~6 



i 2 

1o)z 



(1 - 50 - 3g - jo)z 



(34) 



5.2 Kinematic approach 

A different way of describing the cosmic expansion history is 
by means of the so-called kinematic models that do not use 
quantities from the dynamic description, such as f2 m or w . 
Her e we discuss the models of [Kigar0v fe Multamakil (|2006t ) 
and IGuimaraes. Cunha fc Limal (|2009T ), that are based on 
different parameterisations of the deceleration parameter, q, 
and the jerk parameter, j. The deceleration parameter in 
terms of z is defined as 



la I,, AH(z) 2 )' 



(32) 



where the primed quantity denotes the derivative with re- 
spect to z. Similarly, the jerk parameter is given by 

1 1 d 3 a 



2 ( + > H 2 



(33) 



We consider five realizations here: in the first and sim- 
plest model, Mo, the deceleration parameter is constant, 
q(z) = go. The second model, Mi, is a linear expansion 
of the deceleration parameter q(z) — qo + qiz. Model M2 
has two phases of constant deceleration, separated by an 
abrupt transition redshift, i.e. q(z) = go for z ^ z t and 
q(z) = gi for z > z t . The fourth model, M3 is a constant- 



All these model s make specific predicti ons for the cos- 
mic expansion (see [G uimaraes et al. 1 120091, for the corre- 
sponding expressions of H(z)). They are of particular in- 
terest because no assumption on t he matter-energy c ontent 
present in the Universe is made. iGuimaraes et all HH) 
constrained the different parameters in their models through 
a Bayesi an statistical analysis from the previous Union com- 
pilation (|Kowalski et al.ll200sl l 

The comparison between the expansion history ob- 
tained from Union2 data and the theoretical expansion func- 
tion calculated for each model is shown in Fig. [S] For the 
different parameters we adopt the best-fit values obtained by 
IGuimaraes et all (|2009l l Again, our goal here is to illustrate 
how the method performs and how families of models, with 
completely different parameters, can be confronted with the 
data in order to assess their validity. We note that, for ex- 
ample, model M2 could be rejected based on current SN la 
data because its shape differs greatly from the reconstructed 
expansion history. Also Mo has a different slope which does 
not even fit the data at low redshifts. Models Mi, M3 and 
M4 agree with the SN data within the error bars up to values 
of the scale factor a ~ 0.6 (2 ~ 0.6 6). This is consistent wit h 
the transition redshifts z t given in IGuimaraes et al. I (|2009t ). 
The good agreement may reflect the fact that the models 
were constructed in a way that mimics ACDM at low red- 
shifts. At high redshifts z > z t , however, none of the models 
is consistent with our reconstructed expansion history. 
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6 IMPROVING THE RECONSTRUCTION 
WITH MOCK SN IA DATA 

More data at various redshifts will be obtained in future 
SN la surveys. In order to test for the effects new data may 
have on the reconstructed expansion function H(z), we cre- 
ate a mock data sample of nearby SNe la with fixed Ho- We 
generate a sample of 200 mock SNe la in a ACDM model 
with Q m = 0.274, n d e = 0.726 and H = 70 km/s, with 
an uniform distribution in scale factor from 0.86 to 1. The 
reason for adding 200 SNe is to approximately equal the 
number of nearby and distant SNe contained in Union2. We 
calculate the luminosity distances consistent with equation 
Eq. (O and generat e random errors wit h uncertainties 30% 
smaller than those of lHicken et al.l (|2009T ) . The random num- 
ber is drawn from a normal distribution with mean zero and 
a standard deviation of one. 

The result of the model-independent reconstruction 
based on the Union2 sample including the mock data is 
shown in Fig. [7] The improvement in the reconstruction 
is apparent from the significant decrease in the error bars. 
The reduction of the error at low redshifts was forced by 
constructing the mock data assuming a fixed Ho- This high- 
lights the importance of a well-measured value of Ho for 
cosmological studies. Although SNe la are the best relative 
distance indicators out to redshifts z > 1, they do not pro- 
vide reliable absolute distances. Therefore, the zero point in 
our H(z) reconstruction is not well-determined. Fixing Ho 
from other measurements as a prior in our analysis would 
be beneficial. 

A ACDM model resulting from fl m = 0.274l° °|° 
jAmanullah et al.ll201o]| is shown in Fig. for comparison. 
This model is consistent with our reconstruction within the 
error interval. However, a small difference in slope is ob- 
served, similar to that for the original Union2 sample (cf. 
Fig. [2]). The slight deviation at high redshifts is probably 
due to the additional weight put on the mock data which 
have smaller errors. 

This analysis demonstrates the power of our method in 
testing redshift ranges in which more (and more accurate) 
SN data would help to significantly improve the constraints 
on cosmological models. Such considerations may help with 
the design of future surveys. 

7 TESTING FOR SYSTEMATICS IN SN IA 
DATA COMPILATIONS 

The method of model-independent reconstruction of the ex- 
pansion history of the Universe can in principle be used as 
a tool to study systematics in SN la data. Given, for in- 
stance, two samples from different surveys, a question to 
ask is whether the respective cosmological parameters de- 
rived from both of them predict an expansion history that 
is consistent with the direct data analysis. 

We demonstrate this on the example of the Union2 data 
set which is a collection of data from various SN la sur- 
veys. As noted in Sect. 14.11 there is some tension between 
the reconstructed expansion history and the predictions of a 
ACDM cosmology at intermediate redshifts, although within 
the error bars they are consistent. One possibility is that this 
tension could be due to different characteristics and system- 
atic errors of the sub-samples. In the most extreme case, 



Union2 Sample plus 200 mock nearby SNe la 

1.5 0.7 0.2 0.0 
1801 ' ' ' 1 ' ' ' 1 ' ' ' 1 




60 | i i i | i i i | i i i 

0.4 0.6 0.8 1.0 



Figure 7. The reconstructed expansion rate with 3<x errors ex- 
tracted from the Union2 set, p lus 200 nearby mo ck SNe with 
errors 30% smaller than those of lHicken et al.l 120091) . Our recon- 
struction is compared with a ACDM model assuming the bes t-fit 
cosmological parameters reported by [Amanull ah et al. 

one sub-sample could be responsible for introducing the di- 
vergence from a ACDM model. This can be identified by 
reconstructing the expansion history for the sub-samples in- 
dividually. The two main components of the Union2 data 
set are the ESSENCE and the SNLS samples, as they both 
cover a wide range in redshifts with a relatively large number 
of objects. 

The expansion histories of both samples are compared 
with ACDM models based on the individual best-fit param- 
eters in Fig. \S\ Generally, the error bars for the ESSENCE 
data are wider because they concentrate on a redshift range 
of about 0.6 < a < 0.8 (0.25 < z < 0.66), whereas the 
SNLS data extend up to a ~ 0.5 (z ~ 1). But, apart from 
that, both samples show consistency of the reconstructed ex- 
pansion with a ACDM model. The tension at intermediate 
redshifts is present in both of them. Therefore, the devia- 
tion of the slope of the reconstructed H(z) from the ACDM 
prediction seen in the Union2 data is not introduced by a 
mismatch of the main sub-samples or systematic errors in 
one of them. 



8 CONCLUSIONS 

In this work, we provided further constrains on the ex- 
pansion history of the Universe in a model-independent 
way. The luminosity-distance measurements, obtained from 
SNe la, depend only on space-time geometry, and can be di- 
rectly related to the Hubble function without need of assum- 
ing any dynamical model. To illustrate the power of this ap- 
proach to investigate Dark Energy models, the method was 
used to reconstruct the expans ion history from the Union2 
sample (|Amanullah et al]|2010l ) - the most complete compi- 
lation of SN la distances to date. Furthermore, a comparison 
with ACDM models with different values of the matter den- 
sity, f2 m , and the dark energy density, Qde, was presented. 
For the Union2 sample, we found good agreement between a 
ACDM model with the best fit parameters obtained specif- 
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ESSENCE Sample 



3a reconstruction 

AC DM model: 

Q m = 0.32^5 (best tit tram Wood-Vasey et al. 2007) " 
w = -1 — 




SNLS Sample 



3o reconstruction 

AC DM model: 

Q. m = 0.263±0.042 (best fit from Astier et al. 2006) 
w = -1 




Figure 8. The reconstructed expansion history with 3(7 errors 
for the ESSENCE (upper panel) and SNLS (lower panel) sam- 
ples compared to ACDM models based on the individual best- 
fit parameters. Th e valu e of Q m = 0.32±^ was obtained by 
llWood-Vasev et al.ll2007l) combinin g the ESSENCE, SNLS and 
the gold samples l|Riess et alj|2004l '). In order to get a reliable re- 
construction 57 nearby SN la from literature have been added to 
both samples. 



ically from this sample and the data. However, the slope of 
the curve is not exactly the same and a systematic deviation 
seems to appear at intermediate values of the scale factor. 

General relativity can accommodate the detected ac- 
celeration of the Universe but cannot give a deeper under- 
standing about its cause. Non-standard cosmological models 
have been suggested as alternative explanations for accel- 
eration without dark energy. We made use of the model- 
independent method to quantitatively compare the ability 
of some dark energy models and non-standard cosmologies 
to represent the data. Comparing the reconstruction recov- 
ered from the Union2 sample data with some braneworld, 
f(R), and kinematic models we found that none of the con- 
sidered cosmologies can be rejected rigorously on the basis 
of current SN la data. Particular braneworld models, such 
as the mDGP/DGP with some parameters, however, as well 
as some f(R) and some kinematic models clearly disagree 
with our reconstruction and can be ruled out. Using an opti- 



mal set of functions based on principal component analysis 
instead of an arbitrary basis will strengthen this statement. 

We argue that the model-independent method is also 
a promising tool in checking for biases and discrepancies in 
the SN la cosmology samples. This was tested by perform- 
ing the reconstructions for sub-samples of the Union2 set - 
specifically the ESSENCE and the SNLS data - separately 
and comparing them with the corresponding best-fit ACDM 
models. We find a general agreement between the different 
data and the models, although some tension at intermediate 
redshifts seems to be present in both. 

Although the method as employed here is indicative of 
viability of dark energy models, improvements are neces- 
sary to reach full decisive power. Apart from constructing 
an optimal basis for the series expansion of the luminosity 
distance, the model-dependence induced by the calibration 
of the SNe la as distance indicators has to be taken into 
account . 

An increasing body of high quality data will help to 
reduce systematic uncertainties and will provide a better 
reconstruction of the expansion rate. To demonstrate this, 
we studied the effect that more nearby SNe la will have 
on the reconstruction of H(a). We generated a sample of 
200 SNe at very low redshifts, uniformly distributed in the 
ra nge 0.86 < a < 1 an d with 30% smaller errors than those 
of iHicken et~aH (|2009l ). These mock data were added to the 
Union2 sample and we found a great decrease in the size 
of the error bars. This confirms the fact that not only an 
increased number of distant SNe, but also higher quality 
data and larger number of nearby SNe are indeed needed to 
reconstruct the expansion history more accurately. 

The strength of the method employed here to derive the 
expansion history of the Universe is that it provides a purely 
geometrical test. Contrary to many other ways of analyzing 
cosmological data it does not revert to assumptions on the 
energy contents of the Universe nor its dynamics. It thus 
offers a complementary way of detecting possible system- 
atic effects which could affect the data and be overlooked 
within a traditional analysis based on physically motivated 
parametrization. Moreover, it can be used to discriminate 
between different cosmological models and break the degen- 
eracy in the cosmological parameters. In future work, we 
will consider not only SN la data but also other cosmologi- 
cal probes, such as the angular distances from BAO, in or- 
der to obtain tighter constrains on the expansion history of 
the Universe. Finally, is it worthwhile noting the potential of 
the method for the analysis of possible local inhomogeneities 
through comparison of the expansion history in different di- 
rections. 
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